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Stability  properties  of  an  intense  proton  layer  (P-layer)  immersed 

in  a background  plasma  are  investigated  within  the  framework  of  a hybrid 

model  in  which  the  layer  ions  are  described  by  the  Vlasov  equation,  and 

the  background  plasma  electrons  and  ions  are  described  as  macroscopic, 

cold  fluids.  Moreover,  the  stability  analysis  is  carried  out  for 

frequencies  near  multiples  of  the  mean  rotational  frequency  of  the  layer. 

It  is  assumed  that  the  layer  is  thin,  with  radial  thickness  (2a)  much 

smaller  than  the  mean  radius  (Rq) . Electromagnetic  stability  properties 

are  calculated  for  flute  perturbations  (9/3z=0)  about  a P-layer  with 

rectangular  density  profile,  described  by  the  rigid-rotor  equilibrium 

0 

distribution  function  f^=(m^n^/2ir)  d(U-T)G(vz) , where  n^  and  T are 

constants,  m.  is  the  mass  of  the  layer  ions,  G(v  ) is  the  parallel 
1 z 

velocity  distribution,  and  U is  an  effective  perpendicular  energy 
variable.  Stability  properties  are  investigated  including  the  effects 
of  (a)  the  equilibrium  magnetic  field  depression  produced  by  the  P-layer, 

(b)  transverse  magnetic  perturbations  (<5B^0) , (c)  small  (but  finite) 
transverse  temperature  of  the  layer  ions,  and  (d)  the  dielectric 
properties  of  the  background  plasma.  All  of  these  effects  are  shown  to 
have  an  important  influence  on  stability  behavior.  For  example, 
for  a dense  background  plasma,  the  system  can  be  easily  stabilized  by 
a sufficiently  large  transverse  temperature  of  the  layer  ions. 
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I.  INTRODUCTION 


Field-reversed  ion  layers  and  rings  have  received  considerable 
recent  attention  as  magnetic  confinement  configurations  for  fusion 


plasmas. 


1-8 


Such  layers  and  rings  are  likely  subject  to  various 


9-13 

macro-  and  micro instabilities.  For  example,  recent  theoretical 

11-13 

studies  of  the  negative-mass  stability  properties  of  a weakly 

diamagnetic  ion  layer  embedded  in  a background  plasma  predict  instability 

for  perturbations  with  frequency  near  harmonics  of  the  layer  rotational 

13 

frequency.  These  studies  have  been  carried  out  for  a low- intensity 

2 2 

ion  layer  characterized  by  v<<1,  where  v=N,  e /m.c  is  Budker's  parameter 

rR  „ b 1 

I c 0 

for  the  layer  ions,  and  N^=2irj  dr  r n^(r)  is  the  number  of  ions 
per  unit  axial  lnegth.  A more  general  stability  analysis  is  required 
to  investigate  stability  properties  for  an  intense  field-reversed 
ion  layer  characterized  by  v>>1. 

This  paper  develops  a hybrid  theory  of  the  negative-mass  instability 

for  intense  ion  layers  with  arbitrary  degree  of  field  reversal.  The 

13 

present  work  extends  the  previous  self-consistent  theory  of  the 
negative-mass  instability  developed  for  v<<1.  The  analysis  is  carried 
out  within  the  framework  of  a hybrid  (Vlasov-fluid)  model  in  which 
the  layer  electrons  and  background  plasma  electrons  and  ions  are 

described  as  macroscopic,  cold  fluids  immersed  in  an  axial  magnetic 

\ 

0 

field  Bz(r)^z,  and  the  layer  ions  are  described  by  the  Vlasov  equation. 

We  assume  that  the  layer  is  thin  [Eq.  (1)],  i.e.,  the  radial  thickness 
(2a)  of  the  layer  is  small  in  comparison  with  the  mean  radius  Rq. 
Equilibrium  and  stability  properties  are  calculated  for  the  specific 
choice  of  ion  layer  distribution  function  [Eq.  (2)] 

m.n, 

- -IT  5(u'T)G(vz)  • 
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f£<HA-a>0P0,  vz) 


where  n^,  u>0 , T are  constants,  G(vz)  is  the  parallel  velocity 

distribution,  H,  is  the  perpendicular  energy,  PQ  is  the  canonical 

y 

angular  momentum,  and  U is  the  effective  energy  variable  defined  in 
Eq.  (3). 

One  of  the  important  features  of  the  equilibrium  analysis 
(Sec.  II)  is  that  the  equilibrium  distribution  function  in  Eq.  (2) 
corresponds  to  a sharp-boundary  density  profile  [Eq.  (9)],  with  uniform 
angular  velocity  profile  over  the  layer  cross  section  [Eq.  (12)], 
and  nonzero  transverse  temperature  [Eq.  (13)].  Moreover, 
defining  the  magnetic  compression  ratio  n by  n=Bz(r=R^) /Bq  [Eq.  (16)], 
where  Bz(r=R^)  and  Bq  are  the  axial  magnetic  fields  at  the  inner  and 
outer  surfaces  of  the  layer,  we  find  [Eq.  (26)] 

v=(l-n)/(l+n)  , 


for  a thin  layer  with  a<<R^.  Because  v>0,  it  is  important  to  note  from 
Eq.  (26)  that  the  allowable  range  of  compression  ratio  is  given  by  -l<n<l. 
Moreover,  in  order  to  produce  a sizeable  field  depression,  a relatively 
large  value  of  Budker's  parameter  (v*l)  is  required.  We  also  emphasize 
that  the  rigid-rotor  distribution  function  in  Eq.  (2)  has  an  associated 
spread  in  canonical  angular  momentum  (A)  [Eq.  (32)],  which  plays  an 
important  role  in  determining  stability  behavior.^ 

The  electromagnetic  stability  analysis  in  Secs.  III-V  includes  the 
effects  of  (a)  the  equilibrium  magnetic  field  depression  produced 
by  the  P-layer,  (b)  transverse  magnetic  perturbations  (5B^0) , (c)  small 
(but  finite)  transverse  temperature  of  the  layer  ions,  and  (d)  the 
dielectric  properties  of  the  background  plasma.  The  analysis  is  carried 
out  within  the  framework  of  the  linearized  Vlasov-fluid  and  Maxwell 
equations,  assuming  that  all  perturbed  quantities  are  independent  of 


....  _ - _ „ . . . 


g3  + g2  - n | 


where  G is  the  normalized  Doppler-shifted  eigenf requency  defined  in  Eq. 

(88),  and  the  parameters  C and  A are  defined  in  Eqs.  (87)  and  (91), 
respectively.  In  Eq.  (90),  [Eq.  (87)]  is  an  oscillatory  function 
of  plasma  density.  However,  the  value  of  s',  averaged  over  each  period, 
is  an  increasing  function  of  plasma  density.  We  therefore  conclude 
from  Eq.  (90)  that  the  system  is  completely  stabilized  if  the  plasma  density 
is  sufficiently  high.  The  terms  proportional  to  i,A  in  Eq.  (90)  also  have 
a stabilizing  influence,  thereby  quenching  the  growth  rate  for  sufficiently 
high  Z values.  This  effect  is  most  pronounced  when  the  magnetic 
compression  ratio  q is  close  to  zero.  A similar  stabilization  for  high  Z 
values  has  also  been  demonstrated  for  intense  relativistic  E-layers.^3- 
A numerical  investigation  of  stability  properties  is  carried  out 
in  Sec.  V.C  for  general  values  of  an^  (the  plasma  density  outside  the 
layer).  Several  points  are  noteworthy  in  this  regard.  First,  the 
instability  growth  rate  is  greatly  reduced  as  q-»0 . This  feature  is  evident 
from  Eq.  (90)  for  a=0.  Second,  stability  properties  are  almost  independent 
of  a,  provided  a is  sufficiently  small  (a<0.5,  say).  Finally,  in 
parameter  regimes  where  instability  does  exist,  the  maximum  growth 
rate  can  be  a substantial  fraction  of  ion  cyclotron  frequency 
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II.  EQUILIBRIUM  THEORY 

A.  General  Equilibrium  Properties 

The  equilibrium  configuration  is  illustrated  in  Fig.  1.  It 
consists  of  a space-charge  neutralized  P-layer  (proton  layer)  that  is 
infinite  in  axial  extent  and  immersed  in  a cold,  dense,  background 
plasma.  The  plasma  ions  are  assumed  to  be  singly  charged,  and  the 
mean  radius  and  radial  thickness  of  the  P-layer  are  denoted  by  Rq 
and  2a,  respectively.  The  radius  of  the  cylindrical  conducting  wall 
is  denoted  by  Rc>  The  mean  motion  of  the  P-layer  is  in  the 
azimuthal  direction,  and  the  applied  and  self  magnetic  field  provides 
radial  confinement  of  the  layer  ions.  As  shown  in  Fig.  1,  we  introduce 
a cylindrical  polar  coordinate  system  (r,0,z)  with  z-axis  coinciding 
with  the  axis  of  symmetry;  r is  the  radial  distance  from  the  z-axis, 
and  0 is  the  polar  angle  in  a plane  perpendicular  to  the  z-axis. 

The  following  are  the  main  assumptions  pertaining  to  the  equilibrium 
configuration: 

(a)  Equilibrium  properties  are  azimuthally  symmetric  (3/90=0) 
and  independent  of  z. 

(b)  The  radial  thickness  of  the  P-layer  is  much  smaller  than  its 
major  radius,  i.e.. 


a«RQ  . (1) 

(c)  The  background  plasma  electrons  (j=e)  and  ions  (j=i),  and 

the  layer  electrons  (j*e')  are  treated  as  macroscopic  cold  fluids  (T^=0) 

In  equilibrium,  these  fluids  are  assumed  to  be  stationary  with  zero  net 

axial  motion  [V^.(r)=0]  and  zero  net  azimuthal  motion  [vj? . (r)*w. (r)r=0] 
zj  eJ  J 

for  j=e,i,e' . 
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(d)  The  background  plasma  equilibrium  is  assumed  to  be  electrically 
neutral  with  n?(r)=n^(r).  In  addition,  the  equilibrium  charge  density 
of  the  layer  ions  (j=b)  is  neutralized  by  the  layer  electrons  (j=e') 
with  n^(r)=n^, (r) , and  the  equilibrium  radial  electric  field  is 
equal  to  zero,  E^(r)=0. 

For  the  layer  ions,  any  distribution  function  f^(^c,^)  that  is  a 
function  of  the  single-particle  constants  of  the  motion  in  the 
equilibrium  fields  is  a solution  to  the  steady-state  (3/3t=0)  Vlasov 
equation.  For  present  purposes,  we  consider  the  class  of  rigid-rotor 
Vlasov  equilibria  described  by^ 


*;<HrV9.  V ■ -fir  «<»-wc<v,)  , 


where  n^,  , and  T are  constants,  G(vz)  is  the  parallel  velocity 


distribution  with  normalization  dv  G(v  )=1,  and  the  effective  energy 

J _oo  z z 

variable  U is  defined  by 


B-Hi-»9P0-«1R^+<e/c)Ro»9A"  (R0) 


In  Eq.  (3),  H±  is  the  perpendicular  energy 

Hi=(mi/2) (v^+Vq)  , 


and  PQ  is  the  canonical  angular  momentum 


Vrtmive+(e/c)Vr)]  • 


Here  e and  m^  are  the  proton  charge  and  mass,  respectively,  c is  the 
speed  of  light  in  vacuo,  and  v and  vQ  are  the  radial  and  axial  velocities 


of  a layer  ion.  The  0 -component  of  the  equilibrium  vector  potential, 
Aq(t),  is  to  be  calculated  self-consistently  from  the  steady-state 


— 


Maxwell  equation. 


313  A0,  , 

37  717  rAe(r) 


4Tre  0,  x 0,  , 
— nb(r)V0(r) 


(6) 


where  n^(r)  is  the  local  ion  layer  density. 


n£(r)=  d3v  f°(HA-o.0P0,  v() 


(7) 


and  V (r)  is  the  mean  azimuthal  velocity  of  an  ion  layer  fluid 
0 


element , 


V°(r)-< 


d v v0fb)/( 


0v„f,  3 


d v 9 • 


(8) 


Substituting  Eq.  (2)  into  Eq.  (7),  we  find  that  the  ion  layer 
density  is  given  by  (Fig.  2) 


0/  \ 
nb(r) 


% ’ Rl<r<R2  ’ 


(9) 


0 , otherwise  , 


where  and  R£  are  the  extremes  of  the  interval  on  which  the  inequality 


<|»  (r)>0 


(10) 


is  satisfied.  [That  is,  R^  and  R^  are  determined  from  i|»(R^)='KR2)"0.  ] 


In  Eq.  (10),  the  envelope  function  t|/(r)  is  defined  by 


4)(r)  = (miWg/2)  (r2-RQ)  + (e/c)o)0  [rA0  (r)-RQAg(Ro)  ]+T  . (11) 


From  Eqs . (2)  and  (11),  ^(r)=(m  /2) [v2+(v  -rw  )2]  is  the  (r-0) 

I r t)  d 


kinetic  energy  of  a layer  ion  in  a frame  of  reference  rotating  with 
angular  velocity  ma . In  a similar  manner,  the  azimuthal  velocity 

U 

profile  associated  with  the  equilibrium  distribution  function  in  Eq. 
(2)  can  be  expressed  as 


■ 
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VQ°(r)  = 


:roi„ 


Rl<r<R2 


(12) 


Equation  (12)  corresponds  to  a rigid-rotor  angular  velocity  profile 
over  the  layer  cross  section.  Defining  the  effective  temperature 
T?(r)  for  the  transverse  motion  of  the  layer  ions  by 


n£(r)T°(r)= 


d3v[v^+(v0- 


we  find  from  Eqs.  (2)  and  (11)  that  the  transverse  temperature  profile 
can  be  expressed  as 


T][(r)-*(r) 


R1<r<R2 


(13) 


Evidently,  the  envelope  function  i//(r)>0  defined  in  Eq.  (11)  is 
identical  to  the  transverse  temperature  profile. 

We  note  from  Eq.  (9)  that  the  density  profile  has  sharp  radial 
boundaries  at  and  R2>  Substituting  Eq.  (12)  into  Eq.  (6).  we 
find  that  the  equilibrium  axial  magnetic  field  within  the  P-layer 
(R^<r<R2)  can  be  expressed  as 

B^(r)=Bo+(2ire/c)o)0nb(R2-r2)  , (14) 


where  Bq  is  the  axial  magnetic  field  at  rH^.  The  corresponding 
equilibrium  vector  potential  A^(r)  is  given  by 


rA2(r)-V9(V+(B0/2)(t2-R5) 


+(Tre/2c)m0nb(r2-R2)  [2R2-(r2+R2)  ] 


(15) 


For  notational  convenience  in  the  subsequent  analysis,  we  introduce 
the  magnetic  compression  ratio  n defined  by 


I 
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1 > 


L 


n=Bz(r=V/B0  * 


(16) 


which  characterizes  the  change  in  axial  magnetic  field.  We  further 


introduce  Budker's  parameter  for  the  ayer  ions 


R 

r c 


v=2ir(e  /m^c) 


dr  r 


0 


n^(r) 


(17) 


Substituting  Eq.  (14)  into  Eq.  (16)  and  making  use  of  Eq.  (17) 
gives 


Wq/id  =(n-l)/2v  , 

D Cl 


(18) 


where  u^^eBg/m^c  is  the  ion  cyclotron  frequency  at  r=R2* 


The  axial  magnetic  field  profile  is  illustrated  in  Fig.  3. 

Equations  (11)  and  (15)  completely  determine  the  functional  form 
of  the  envelope  function  ’Kr) . We  note  from  Eq.  (15)  that  equilibrium 
layer  density  (the  term  proportional  to  n^)  can  generally  have  a large 


nonlinear  influence  on  the  location  of  the  radial  boundaries  R^  and 


R£.  A sketch  of  ip(r)  versus  r is  illustrated  in  Fig.  4.  Thus  far, 


Rq  has  been  introduced  in  the  analysis  as  an  unspecified  constant 


parameter  in  Eq.  (3).  Without  loss  of  generality,  we  choose  R 


0 


to  correspond  to  that  radius  where  ^(r)  passes  through  a maximum 


in  the  interval  R^<r<R2  (Fig.  4), 


£♦(*) 


=0  . 


(19) 


r=Rr 


Substituting  Eq.  (11)  into  Eq.  (19)  we  find 


a)0+(e/mic)Bz(Ro)=O  , 


(20) 


where  use  has  been  made  of  Eqs.  (14)  and  (15). 


i 
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Equation  (20)  is  simply  a statement  of  radial  force  balance  on  an 
ion  layer  fluid  element  at  r=Rg.  Further  general  equilibrium  properties 
associated  with  the  distribution  function  in  Eq.  (2)  are  discussed  in 
Ref . 4. 


we  readily  determine  R^Rg-a  and  R2=Rg+a  from  iKr)=0  and  Eq.  (21). 
It  is  also  noteworthy  from  Eq.  (22)  that  the  betatron  frequency 

is  directly  proportional  to  u^Bq  for  the  rigid-rotor  equilibrium 
described  by  Eq.  (2).  This  is  considerably  different  from  the 
result  obtained  in  Ref.  12  for  the  choice  of  distribution  function 
in  which  all  of  the  layer  ions  have  the  same  values  of  canonical 
angular  momentum  and  the  same  value  of  energy. 


For  a thin  P-layer,  it  is  straightforward  to  show 


12 


from  Eq.  (21).  Substituting  Eq.  (11)  into  Eq.  (24)  and  making  use 
of  definition  in  Eq.  (16)  gives 

V“ci=-(1+n)/2  , (25) 

which  determines  the  rotational  frequency  ujq  in  terms  of  the  compression 
ratio  q.  Moreover,  substituting  Eq.  (25)  into  Eq.  (18),  we  find 

v=(l-q) / (1+q)  . (26) 

Equation  (26)  is  one  of  the  most  important  equilibrium  results 
pertaining  to  the  choice  of  distribution  function  in  Eq.  (2) , and 
several  points  are  noteworthy.  First,  Budker’s  parameter  v for  the 
layer  ions  can  be  expressed  exclusively  in  terms  of  the  magnetic 
field  compression  ratio  q.  Second,  since  v>0,  the  allowed  range  of 
q is  given  by 


-l<q<l  . (27) 

Moreover,  Budker's  parameter  rapidly  increases  to  infinity  as  the 
compression  ratio  approaches  minus  unity  (q-*~  1) . In  this  regard, 
strictly  speaking,  it  is  not  possible  to  achieve  complete  field 
reversal  (q=-l)  for  the  choice  of  distribution  function  in  Eq.  (2). 

We  note  that  v=l  corresponds  to  zero  compression  ratio  (q=0) . Therefore, 
in  order  to  produce  a significant  field  depression,  we  conclude 
that  a layer  with  a reasonably  high  value  of  Budker’s  parameter  (v>l) 
is  necessary. 

Of  considerable  interest  for  experimental  application  is  the 
relationship  between  the  applied  axial  magnetic  field  B before  injection 
of  the  layer,  and  the  external  magnetic  field  Bq  after  injection.  In 
this  regard,  we  assume  that  the  magnetic  flux  inside  the  cylindrical 


13 


conductor  is  conserved,  i.e., 


R 

r c 


2ir 


rdr  B^(r)=irR2B  . 
z c 


(28) 


For  a thin  layer,  the  axial  magnetic  field  profile  can  be  approximated 
by 


B^(r)  = Bq 


0<r<RQ  , 


1 , Ro<r<Rc  , 


and  it  is  straightforward  to  show  from  Eq.  (28)  that 


B/B0=l-(l-n)(R0/Rc)  • 


(29) 


For  B/Bq>0,  we  conclude  from  Eq.  (29)  that  the  magnetic  compression  ratio 
is  restricted  to 


n<(Rc/RQ)  -1  . 


(30) 


Finally,  we  conclude  this  section  by  noting  that  the  choice  of 
equilibrium  distribution  function  in  Eq.  (2)  yields  a spread  in 
canonical  angular  momentum  P_.  Defining  the  average  canonical 
angular  momentum  Pq  as 

po=miRoV(e/c)RoAe(Ro)  * 


it  is  straightforward  to  show  from  Eqs.  (2)  and  (11)  that 

| PQ-P0 | =miR0 [2T/m.-vr(R0) 2 ] 1/2 , (31) 

where  v (R  ) is  the  radial  velocity  of  a layer  ion  at  r*Rn.  For  a 
r (J  u 

given  transverse  temperature  T,  those  particles  with  vr(Rg)=0  have 
the  maximum  deviation  in  canonical  angular  momentum  from  the  mean  value 
Pg.  Therefore,  the  maximum  canonical  angular  momentum  spread  A can  be 
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III.  PERTURBED  CHARGE  DENSITY  FOR  A THIN  ION  LAYER 

In  the  present  article,  we  consider  perturbations  with  very  long 
2 2 

axial  wavelength  (kzRg<<l) , and  approximate 

k -0  (33) 

z 

in  the  subsequent  analysis.  (Here  k^  is  the  axial  wavenumber.)  It 
is  assumed  that  all  perturbed  quantities  6i|i(^c,t)  can  be  Fourier- 
decomposed  according  to 


6i//(£,t)=i*4(r)exp{i(£.e-(i)t)},  lmw>0  , (34) 

where  uj  is  the  complex  eigenfrequency,  and  l is  the  azimuthal  harmonic 
number.  Moreover,  the  stability  analysis  is  carried  out  for  eigen- 
frequency near  multiples  of  the  mean  P-layer  rotational  frequency,  i.e., 

0)-£Cl)g  . 

The  Maxwell  equations  for  the  perturbed  fields  become 

V*|(£)=i(u/c)B(£)  , 

Vx^(x)  = (41r/c)^(^)-i(a)/c)E(x)  . 

Within  the  context  of  Eq.  (33),  the  perturbed  electromagnetic  fields 
can  be  expressed  as 


&t<r)'5M(r)«r+ite(r)«e  • 

«*<r)‘Su(r>i  • 


(35) 


Fourier  decomposing  the  perturbed  fields  according  to  Eqs.  (34)  and 
(35),  the  and  V*E  Maxwell  equations  can  be  expressed  as 


* ? ? ? ? ? -1 

E*r(r)“£  (“  r /c  } <3/3rH 

. .2,222  2,-1-  , , 

-4irir  u(u)  r -l  c ) J (r)  , 

£r 

B£z(r)=-Jla)r(to2r2/c-Jl2c)_1(3/3r)i 

+4irirJl(u)2r2/c-£2c)  2J  (r)  , 

ir 


where  the  function  $(r)  is  defined  by 


<j>(r)  = irE£0(r)/£  , 


and  (r)  is  the  perturbed  radial  current  density. 

After  some  straightforward  algebra,  the  perturbed  ion  layer 
distribution  function  can  be  expressed  as 22 

8fb  r° 

fbi(r’&=e  ITu"  dx  exP^-i[wT-£(0'-e) ]} 

—oo 

r , , , , 

x (9--V  (iU  + — ~ (EZr  + ^ , 

where  x=t'-t.  In  obtaining  Eq.  (38),  use  has  been  made  of  3U/3v= 

% 

m^[vrer+(v0-rtD0)ee] , where  e^  and  eQ  are  unit  vectors  in  the  r and  6 
directions,  respectively.  In  Eq.  (38),  the  trajectories,  x'(t') 
and  v'(t')  satisfy 

'V* 

’ JFZ'  m^cX'*h>W  » 

where  x'(t'=t)=£,  and  v'(t'=t)=v.  The  term  E^+VgB^/c  in  Eq.  (38) 
can  be  further  simplified.  Making  use  of  Eq.  (36)  we  obtain 


+ tc  vrt- -h  * • 


Since  the  eigenf requency  of  the  perturbation  is  given  approximately 
by  ii)=£<D0,  Eq.  (39)  can  be  approximated  by 


Eir^+v0BS,z^r^c*_^3/,9r^r^ 


The  evaluation  of  the  orbit  integral  in  Eq.  (38)  is  generally 
very  complicated.  However,  for  present  purposes,  we  assume  low- 
frequency,  long-wavelength  perturbations  characterized  by 


l 1 2 2 

|u>HU)q|  «u>r  , 


&a/Ro<<a)r/|u)0|  , 


where  a is  the  half-thickness  of  the  layer  defined  in  Eq.  (23). 
Within  the  context  of  Eq.  (41),  it  is  valid  for  a thin  layer  to 
approximate 


(j>(r')=(j)(R0)  , 


9 ,=0+[toQ+(v1/RQ)cosa]T  , 


where  vA  is  the  perpendicular  speed  in  a frame  of  reference  rotating 
with  angular  velocity  w0.  Moreover,  vA  is  related  to  the  effective 
energy  variable  U by 

U=(mi/2)[v^(r-RQ)2]  , (4 

with  v0-rw0=vAcosct  [see  Eqs.  (3),  (11),  and  (21)].  Furthermore, 
the  perturbed  ion  layer  distribution  function  can  be  simplified  as 
eVj.  0 

fb«.^r,^)*1A  R — cosct  3U-  dT<J>exp{-i[a)T-£(9 '-0)  ] } . (4 

0 — 00 


Substituting  Eq.  (42)  into  Eq.  (44)  gives 


fb  A £(vA/R0)cosa 


fb£(r’^"  6 3U  4>(R0)  (w-£a>e)-Uvx/R0)cosci  * 
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IV.  ELECTROMAGNETIC  STABILITY  PROPERTIES 


A.  General  Eigenvalue  Equation 

As  discussed  at  the  beginning  of  Sec.  II,  the  background  plasma 
components  (j=e,i)  and  layer  electrons  (j=e')  are  treated  as  cold  (Tj-K)) , 
macroscopic  fluids  immersed  in  an  axial  magnetic  field  B^(r)^  , 
where  B^(r)  is  approximated  by 


B"(r)-Bo 


0<r<RQ  , 


1 ’ R0<r<Rc  ’ 


for  a thin  layer.  The  momentum  transfer  equation  and  the  continuity 
equation  for  each  cold-fluid  component  can  be  expressed  as 

(!f  + -^(«  + ^r)  • 


(3/3t)nj+X.(njVj)-0  , 

where  n^C^t)  is  the  density,  Vj(^c,t)  is  the  mean  velocity,  and  e^ 
and  in.  are  charge  and  mass,  respectively,  of  a particle  of  species  j. 

As  illustrated  in  Fig.  5,  the  stability  properties  of  the  layer- 
plasma  system  are  investigated  for  perturbations  about  the  step-function 
plasma  density  profiles  specified  by 


nj(r,’”P 


1 , 0<r<RQ  , 


a , R0<r<Rc  , 


for  j*e,i.  Here  a is  an  arbitrary  constant.  To  make  the  stability 
analysis  tractable,  the  following  simplifying  assumptions  are  also  made. 


I 


20 
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(a)  The  P-layer  is  immersed  in  a dense  plasma  background,  i.e., 


Vnr 


(51) 


which  is  easily  attainable  in  the  parameter  regimes  of  experimental 
interest. 

(b)  The  perturbed  charge  density  of  the  layer  ions  can  be  represented 
by  the  surface  charge  contribution  given  in  Eq.  (47),  which  is 


valid  for  a sufficiently  thin  layer  (a<<R^)  and  dense  background 


plasma  (n^<n^) . 

(c)  The  influence  of  current  and  charge  perturbations  associated 
with  the  layer  electrons  are  neglected  in  the  stability  analysis. 
Although  this  is  a good  approximation  for  an  infinitesimally  thin  layer 
with  a/R^-^O,  we  expect  some  modifications  to  the  stability  behavior 
associated  with  finite  layer  thickness. 


For  perturbation  with  kz=0,  Eq.  (49)  can  be  linearized  to  give 


-itonj(r)  + ^ [rn*?(r)Vjr(r)  ]=-i  n° (r)V.Q(r)  , 


-i«Vjr(r)-ejUcj(r)Vje(r)  = ^ Vr)  , 


(52) 


-lttVj0(r)+ej«cj(r)Vjr(r)-i  ^ £ *(r)  , 


where  1 is  the  azimuthal  harmonic  number,  j*e,i  denotes  plasma  species, 


and  use  has  been  made  of  . (r)=V^ . (r)*0  for  j-e,i.  In  Eq.  (52), 

zj  0J 


£j*sgn  e^  , 


“cj(r)"“cj  (“eB0/mjc)x 


n , 0<r<RQ  , 


(53) 


1 * R0<r<Rc  * 


is  the  cyclotron  frequency,  V. (r)  and  n. (r)  are  the  perturbed  fluid 

''O  j 


J.  . . V . 


J 
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velocity  and  density,  and  the  abbreviated  notation  x^(r)  = i^^(r)  has  been 
introduced  for  the  perturbation  amplitudes.  The  perturbed  radial 
electric  field  E^Cr)  in  Eq.  (36)  is  calculated  self-consistently 
from  Eq.  (52)  and  the  definition  of  radial  current  density 


J_(r)=  I e n?(r)V  (r)  . 
j=e,i  2 2 jr 


Defining 


we  obtain 


5=l-(wr/J!.c)  , 


2 to  .(r) 

w (if)  I -V 


x — E (r) 
2 r 


. li  _ 5 V)(r)  £j  ~cJ(r)  A 

^ 3r  * \ lc'  r 2 0)  r 

3 Vi 


2 2 2 2 

where  Vj=to  (r)  and  u»  (r)  is  the  plasma  frequency-squared  defined 


2 * 2 2 
to  . (r)=to  . (=4re  n /m.)x 
PI  PJ  P J 


1 , 0<r<RQ  , 


a , R <r<R 
0 c 


Poisson's  equation  for  the  perturbed  electric  field  can  be 


expressed  as 


7 [rEr(r)]  + ^ <Kr)-4ira  6(r-RQ)+4irj[e  n.(r)  , 
r j 2 2 

where  the  surface  charge  density  is  defined  in  Eq.  (47). 
Eliminating  V^r(r)  and  V^g(r)  from  Eq.  (52)  in  favor  of  n^ (r) , 
and  substituting  n^ (r)  into  Eq.  (57)  gives 

i a ( f - ) o2  _ (rh  - 

7 37  r1'  Vr>  +~  -"V  »(r> 

l ^ i v ; ) r K j V. 


1»  «!<*> 
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r l^r  r Jj  ^ ( 


u)2.  (r) 


y-  ui“cj<r)!i 


2 2"  2 
2£  c v*/Rg 

2 2 2 2 ^ * 

Mwe)  -JlV(a/R0)Z  0 


where  v=4ire  n^R^a/m^c  is  Budker's  parameter  for  the  layer  ions.  It  is 
useful  to  introduce  the  abbreviated  notation 


w2 . (r) 


S1(r,ai)=^  1 , 

j Vi 

2 

<o  . (r)  e.w  .(r) 

s2<r'“)*?  -nT~  — 

vi 

J UJ 


Eliminating  E^Cr)  from  Eqs.  (55)  and  (58),  and  making  use  of  the  definitions 
in  Eq.  (59)  gives  the  eigenvalue  equation 


I i-4  !l  r ii)  fill 

r 9r\4  S3  9r/  r 4 1 S3  3r 
2 q ( 2 

(sf-si+i)]+ 

!c  , 

2 I ur/wr  A X/  r,  \ 

sj  l1+5te)  J |+4”«b«(r-R0) 


where  is  defined  in  Eq.  (47). 

The  eigenvalue  equation  (60)  is  fully  electromagnetic  and  has 

been  derived  with  no  a priori  assumption  on  the  relative  strengths  of 

the  transverse  magnetic  and  electric  perturbations.  If  we  formally 

2 2 2 2 

take  the  electrostatic  limit  in  Eq.  (60)  with  r w /Si  c -K),  then  £->-1 

[Eq.  (54)],  S3-l  [Eq.  (59)],  Er--3$/9r  [Eq.  (55)],  and  Eq.  (60) 
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simplifies  to  give  the  familiar  form 
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7 h (sir  If)  - h V 


(61) 


7 * h S2+4Trab6(r-RQ) 


Of  course,  strictly  speaking,  the  electrostatic  eigenvalue  equation  in 
Eq.  (61)  is  valid  only  for  a low-beta  layer-plasma  system  with  negligible 
magnetic  field  depression. 


B.  Approximate  Eigenvalue  Equation 
For  a layer-plasma  system  with  arbitrary  degree  of  field  reversal, 
the  electrostatic  eigenvalue  equation  (61)  is  not  valid  and  it  is 


necessary  to  make  use  of  Eq.  (60).  For  the  low-frequency  (oo^JZ-to. ) 

0 


perturbations  considered  here,  it  follows  that 

i 2 


rnr 


£c 


_ 2 2 

Ro“e 

'V'  — — <<  1 . 

a.  2 * 


(62) 


and  hence  that  E,  can  be  approximated  by 


5=1 


(63) 


in  Eq.  (60).  On  the  other  hand,  for  arbitrary  degree  of  field 


2 2 2 

reversal,  it  is  necessary  to  retain  terms  proportional  to  (lor/Hc) /v^  • 


We  therefore  approximate  Eq.  (60)  by 
-S,  „ 2 S 


15_A  3£\  l*  al  . /u>r\2/S2  _ \| 

r 3r(s  r 3r/‘  2 S 1 " ( He  ) VS  " Sl)l 

'3  r 3 1 ' J 

" 7 ♦ ’ 


(64) 


where  and  S£  are  defined  in  Eq.  (59),  and 


2 a)  (r) 

*,  - 1 + (!l)  * -V 

j V. 


j 


J 


V — 


(65) 


The  eigenvalue  equation  (64)  is  generally  difficult  to  solve 
analytically.  However,  a formal  dispersion  relation  that  determines 
the  complex  eigenfrequency  w can  be  obtained  in  a relatively 
straightforward  manner.  Since  the  perturbed  azimuthal  electric 
field  EiQ(r)  is  continuous  across  the  layer  (r=R^) , the  function 
<Kr)  Is  also  continuous  at  r=Rg.  A further  boundary  condition  on 

A A 

4>(r)  in  Eq.  (65)  is  determined  from  the  discontinuity  of  (3<j>/3r)  at 
t^Rq.  For  convenience  of  the  subsequent  analysis,  we  define  the  wave 
admittance  b,  as 


b_=(r3<t./3r)r=R-/i$(R0)  , 


b+»-(r3<t>/3r)r=R+/£<f>(R0)  , 

+ 

where  RQ  denotes  lim5_^>+(R0±6) . Multiplying  Eq.  (64)  by  r and 
integrating  with  respect  to  r from  R^-d  to  RQ+d  (with  d-*0+)  gives 

2^c2v/R 2 

2 — 2~2 2 = “D  * < 

(u»-lwQ)  -4  u^a/Rg)^ 

where  D(to)  is  defined  by 

„-i,  x b+s1(Ro,<il)+s2(Ro,“) 

D (to)  = — 

S3(R0,to) 

b_S^  (RQ  ,m)  "“S2  ( Rq  * to) 

+ _ , 

S3  (Rq  » to) 

and  use  has  been  made  of  the  definition  of  [Eq.  (47)]. 

Equation  (67)  constitutes  the  desired  dispersion  relation 
for  the  complex  eigenfrequency  w.  In  order  to  evaluate  closed 
expressions  for  b+,  however,  we  emphasize  that  it  is  necessary  to 
solve  Eq.  (64)  for  the  eigenfunction  $(r).  Although  this  generally 
requires  a numerical  analysis  of  Eq.  (64),  in  Sec.  V analytic 


w 


solutions  for  <j>(r)  are  obtained  in  two  limiting  regimes  of  experimental 
interest. 


In  concluding  this  section,  it  is  instructive  to  solve  Eq.  (67) 
iteratively  for  eigenfrequency  go  in  the  vicinity  of  £ooQ.  If  it  is  valid 

to  approximate  go-£w  on  the  right-hand  side  of  Eq.  (67),  we  obtain 

0 

2 2 

(w-i lw0)2=2V  \ - 21  ~~  DUu>0)  . (69) 

R0  R0 

2 2 2-',  2 

Making  use  of  wr=v^/a  “21/m^a  [Eq.  (23)],  it  follows  from  Eq.  (69) 
that 


2i.vD(ftu).)  >v2/ c2 
0 1 


(70) 


is  a necessary  and  sufficient  condition  for  instability. 

Equation  (69)  constitutes  one  of  the  main  results  of  this  paper 
and  can  be  used  to  investigate  stability  properties  for  a broad  range 
of  system  parameters.  In  this  regard,  we  emphasize  that  Eq.  (69) 
has  been  derived  with  no  a priori  assumption  regarding  the  size  of  the 


electromagnetic  coupling  parameter 

2 2 


K 


-(£)? 


(71) 


For  general  value  of  K,  the  eigenvalue  equation  (64)  must  be  solved 

numerically  to  determine  the  wave  admittances  b+  [Eq.  (66)].  However, 

in  the  limiting  regimes  where  | K | < < 1 or  | K | >> 1 , closed  analytic  expressions 

for  <t>(r)  (and  hence  b+)  can  be  obtained  in  a straightforward  manner. 

For  a low-density  background  plasma  consistent  with  | K | < < 1 , the 

electrostatic  approximation  is  valid  and  the  corresponding  stability 
13 

properties  have  been  investigated  previously  by  the  authors  for  a 
weakly  diamagnetic  configuration.  In  the  subsequent  analysis,  we  therefore 


investigate  stability  properties  for  a high-density  background  plasma 
with  |k|»1  for  0<r<R_,  and  arbitrary  degree  of  field  reversal. 
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V.  STABILITY  ANALYSIS  FOR  A HIGH-DENSITY  BACKGROUND  PLASMA 

A.  Eigenvalue  Equation 


In  this  section,  we  obtain  closed  expressions  for  the  wave  admittance 
b+  [Eq.  (66)],  and  the  results  are  used  to  investigate  the  dispersion 
relation  in  Eq.  (67).  Making  use  of  5=1  and  Eqs.  (53)  and  (56),  it  is 
straightforward  to  show  from  Eq.  (64)  that  the  eigenfunction  <J>(r) 
satisfies 


1 

3 3 

2 / 

1 . ±\l_ 

*2 

CM 

cn  ; 
csl 

r 

1 Ml 

r 

3r  r 3r 

r { 

S3 ) 3r 

2 

r 

l \ 

s3y  j 

+ (t)2  (17  * si)  • 


at  all  radial  points  except  r=Rg.  In  the  limiting  regimes  where 
|k|«1  or  |k|»1,  the  eigenvalue  equation  (72)  can  be  simplified  to  give 

(r  fr  r fr  ' ^l)  ^(r^=0  • |k|«1  , 

2 (73) 

|K|»1, 

where  use  has  been  made  of  S^l  for  |k|«1,  and  S^>>1  for  |k|»1. 

2 2 

For  an<^  m^/me>>l,  the  quantities  S^  and  in  Eq.  (73) 

can  be  approximated  by 

2 /-  ^ 
to  . (r) 

S2  - „cl(r)  S1  - .cl(t)  v2(r)  • <»> 

2 2 

provided  (1-v)  >>(£ra  /m.)  . [This  is  a very  weak  limitation  on  the 

e i 

range  of  v for  which  the  subsequent  stability  analysis  is  valid. 

It  essentially  allows  for  all  values  in  the  range  0<v<°°,  except 
v*l.]  Substituting  Eq.  (74)  into  Eq.  (73),  we  obtain 


’ 

3 13 

2a)  \ , w2  “pi 

r 3r  r 3r  ' 

r2  ( 

Jlco  ,(r ) 2 2 ( \ 

ci  ' c(o.  (r) 

ci 

<(>(r)=0 
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(75) 


for  I K | »1. 

From  Eq.  (37),  it  is  straightforward  to  express  Eq.  (75)  as 

2 


(r  !r  r fr  ~ “2  + p2)  ®0(r)*°  • 


(76) 


where 


““ pi(r) 

P=-  7l/;Cr ■)’ 


q - U+a 


ci 

2 


1 - 


2io 


£“ci(r) 


1/2 


(77) 


In  obtaining  Eq.  (76),  use  has  been  made  of  the  property  that  <o  . (r) 

P1 

and  to  (r)  are  uniform  except  at  r=R.  [see  Eqs.  (53)  and  (56)].  The 
cl  0 

solution  to  Eq.  (76)  is  a linear  combination  of  J^(pr)  and 
N^(pr),  where  J^(pr)  and  N^(pr)  are  the  Bessel  functions  of  the 
first  and  second  kind,  respectively.  In  the  subsequent  analysis,  it  is 
convenient  to  express  Budker's  parameter  (v)  and  the  angular  velocity 


(go-. ) in  terms  of  the  compression  ratio  q [see  Eqs.  (25)  and  (26)] 
U 


v=(l-n) /(l+n)  , 

to0/wci=-(l+n)/2  . 

B.  Rectangular  Plasma  Density  Profile  (a=0) 


(78) 


As  a simple  limiting  case,  we  examine  stability  properties  in 
circumstances  where  the  plasma  density  outside  the  layer  is  identically 
zero,  i.e.. 
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a = 0 . 


In  the  vacuum  region  (Rg<r<lO  , it  follows  that  ] K ] <<1  and  the  solution 


for  <J>(r)  has  the  simple  form  [Eq.  (73)] 


<f>(r)=»<KR0) 


<R0/r)£-(R0r/Rc^ 


i-ayy 


’ R0<r<Rc 


The  wave  admittance  b+  can  be  determined  by  substituting  Eq.  (80) 


into  Eq.  (66)  which  gives 


V(Rf+Ro*)/(Rc*-Ro*>  • 


For  a dense  background  plasma  core  with  | K | >>1  in  the  range  0<r<Rg, 


the  wave  admittance  b_  is  determined  by  solving  Eq.  (76)  for  <J(r) 


Approximating  in  Eq.  (77)  and  making  use  of  Eqs.  (53)  and  (78), 


we  find  that  the  index  q can  be  expressed  as  (for  0<r<Rg) 


q=[l+Jl2(2+l/n)]1/2 


For  analytic  simplicity,  q is  assumed  to  be  real  in  the  subsequent 


analysis , which  restricts  n to  the  range 


0<n<l  , for  n>0  , 


-Kn<-«.2/(l+2«,2) , for  n<0 


Within  the  context  of  Eq.  (82),  the  physically  acceptable  solution  to 


Eq.  (76)  can  be  expressed  as 


Efl(r)«AJ  (pr),  0<r<R(1  , 


where  A is  an  arbitrary  constant,  p=-(w/c)  (co^/  | n | “ci),  and 
2 1/2 

(Airn^e  /m^)  is  the  ion  plasma  frequency.  Substituting  Eqs.  (37) 


I 
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I 
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1 
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The  quantity  u^a  /Rq=2T/iii^Rq  occurring  in  Eq.  (87)  can  be  expressed  in 
terms  of  the  reversal  parameter  n by  making  use  of  Eqs.  (22),  (23), 
and  (78).  After  some  straightforward  algebra,  we  obtain 


a)V/RQ=^i(a/R0)(l-n2)/4  . (89) 

Substituting  Eqs.  ( 87 ) — (89)  into  Eq.  (86),  and  Taylor-expanding 

(J'/J  ) in  Eq.  (86)  about  oj=Jla)Q,  we  obtain  the  dispersion  relation 
q q o 


+ n 


u 


V_ 

c 


Q + 2\ 


' 2n 

?d+n) 


K." 

2J 


0 , 


(90) 


where 


la  1-n 

Ro  4 


(91) 


The  parameter  £ in  Eq.  (90)  is  an  oscillatory  function  of  plasma 

density.  It  is  evident  from  Eq.  (90)  that  the  parameter  £ plays  a 

decisive  role  in  determining  stability  behavior.  Moreover, 

the  system  can  be  stabilized  by  increasing  the  azimuthal  harmonic 

number  l,  since  the  terms  proportional  to  l in  Eq.  (90)  have  a stabilizing 

influence.  The  stabilization  associated  with  sufficiently  large  l is 

provided  by  the  effective  transverse  temperature  (T)  of  the  layer  ions, 

~ 2 1/2 

and  hence  is  associated  with  the  finite  layer  thickness  a*(2T/m^mr) 

[Eq.  (23)].  Moreover,  shown  in  Eq.  (32),  the  quantity  T is  directly 
related  to  the  canonical  angular  momentum  spread  (A)  for  the  class 
of  rigid-rotor  Vlasov  equilibria  described  by  Eq.  (2).  Therefore,  the 
stabilization  provided  by  the  transverse  temperature  T is  also  associated 
with  the  finite  spread  in  canonical  angular  momentum. ^ This  effect 
is  most  pronounced  when  the  reversal  parameter  n is  close  to  zero  [Eq.  (90)]. 

Specific  stability  properties,  determined  numerically  from  Eq.  (90), 


I 


i 


will  be  discussed  in  Sec.  V.C. 


1 


case 


C.  Stability  Analysis  for  Arbitrary  g 
In  this  section,  we  investigate  stability  properties  for  the 
where  the  plasma  density  outside  the  ion  layer  (RQ<r<Rc)  is  sufficiently 
high  that 


“ R0 


. 2 2 
Z c 


^^Ir^R- 

“ci  0 


=|Kl  D+  » 1 . 


(92) 


r=Rij 

Here  aici=eB0/mic  is  the  ion  cyclotron  frequency,  an  is  the  plasma 

2 2 

density  outside  the  layer,  and  cui  ,=4ran  e /m.  is  the  associated 

Pi  P i 

ion  plasma  frequency-squared.  The  required  solution  to  Eq.  (76)  in 


the  range  RQ<r<Rc  is  given  by 


Ee(r)=B[Jq(pr)Nq(PRc)-Jq(pRc)Nq(Pr) 1 » 


(93) 


where  B is  a constant  and  use  has  been  made  of  the  boundary  condition 


Ee<r=R>°-  Substituting  Eqs.  (37)  and  (93)  into  Eq.  (66)  yields 


b+(w) 


i . 1*0 

i « VpVVpEc)-VpRc)V',V  ' 


(94) 


The  dispersion  relation  can  be  derived  by  substituting  Eqs.  (84) 
and  (94)  into  Eq.  (68).  After  some  straightforward  algebraic  manipulation, 
we  obtain  the  approximate  dispersion  relation 


2oj2v 


(a)-£u0)2-A2a)2(a/Ro) 


2 = ~Ub_+b+)  , 


(95) 


where  the  sum  of  the  wave  admittances,  (b  +b  ) , can  be  approximated  by 


p_R  J’  (p  R ) 

b_(oj)+b+(to)  = —■  J 

q-  - 0 


(96) 


/p+*q\ 

1 * > V<P+V  V^V’VVV  V’+V  ' 

In  Eq.  (96),  p_=-(co/c)  (wpi/ 1 n | ioci)  , p+=-(oi/c)  (a1//2upi/a)ci)«  | n | a1^2p_, 
2 1/2  2 

q_-[l+4  (2+1/n)  ] , and  q+=[l+J,  (2+p)].  In  deriving  Eq.  (95),  use 

has  been  made  of  Eqs.  (65),  (74),  and  (92).  Moreover,  in  obtaining 
the  expressions  for  q+,  use  has  been  made  of  Eq.  (77),  and  ui~^d= 
-£(l+n)toc^/2.  Since  the  eigenf requency  to  is  very  close  to  4co0  , we 
Taylor-expand  the  right-hand  side  of  Eq.  (95)  about  to=4ioQ  , retaining 
terms  to  first  order  in  (u)-4co0) . After  some  straightforward  algebra, 
Eq.  (98)  can  be  approximated  by 
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where  to0  and  ^(a/Rp)  are  defined  in  Eqs.  (78)  and  (89),  respectively. 

The  growth  rate  y=Imto  and  real  oscillation  frequency  Reto  have 

been  obtained  numerically  from  the  cubic  dispersion  relation  (97) 

for  a broad  range  of  system  parameters  n,  to  .R n/c,  a,  4,  and  R„/R  . In 

pi  U 0 c 

the  remainder  of  this  section,  we  summarize  several  features  of  the 
stability  properties  determined  from  Eq.  (97).  From  Eqs.  (41)  and  (78), 
the  real  frequency  Rem  can  be  approximated  by 


Reto=4oje=-4(l+ri)  toci/2  . 


(98) 


The  real  frequency  Reto  determined  numerically  from  Eq.  (97)  is 
very  close  to  the  value  in  Eq.  (98).  In  this  context,  we  only 
present  numerical  results  for  the  instability  growth  rate  y. 
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In  order  to  illustrate  the  dependence  of  stability  properties  on 
the  degree  of  field  reversal,  we  calculate  the  wave  admittances  and 
the  instability  growth  rate.  Shown  in  Fig.  6(a)  is  a plot  of 
(b_+b+)^_£u  versus  n,  for  lm2,  u^^Rq/c-4  and  a«*0.  [In  the  limiting 
case  of  no  plasma  outside  the  ion  layer  (a*0) , the  dispersion  relation 
(97)  reduces  identically  to  the  dispersion  relation  (90)  obtained  in 
Sec.  V.B.]  It  is  evident  from  Fig.  6(a)  that  the  curves  representing 
(b  +b.)  „ become  more  vertical  and  the  distance  between  curves 

- + U)  = X,0)g 

decreases  rapidly  as  n decreases  (or,  equivalently,  as  Budker's 
parameter  v increases).  For  q<0.3,  it  is  readily  shown  that  (b  +b,) 

— + (l)*  J 

can  be  approximated  by 


(b  +b  ) . =(b  ) = 

- + W=i6U)g  - Cl)=£.(i)g 


|4(l-n)  ir  ..  . « 4/2 

tan  ^ — (1  + 4/n) 


(99) 


for  the  parameters  in  Fig.  6(a).  In  obtaining  Eq.  (99),  use  has  been 

1/2 

made  of  Eq.  (77)  and  Jn(x)  = (2/irx)  cos(x-nir/2— ir/4)  for  large  x. 

Equation  (99)  provides  a good  description  of  the  behavior  in  Fig.  6(a). 

In  this  regard,  we  do  not  plot  (b_+b+)  for  n<0.2.  Figure  6(b) 

shows  a plot  of  the  normalized  growth  rate  y/u  ^ versus  q obtained  from 

Eq.  (97)  for  a=0,  Rq/Rc=0.5,  a/RQ=0.05,  and  parameters  otherwise 

identical  to  Fig.  6(a).  Several  points  are  noteworthy  in  Fig.  6(b). 

First,  the  system  is  stabilized  when  the  reversal  parameter  q approaches 

zero.  This  feature  has  been  predicted  analytically  in  Sec.  V.B. 

Second,  the  growth  rate  curve  exhibits  a repetitive  behavior,  with  the 

maximum  growth  rate  for  each  unstable  zone  decreasing  as  q decreases. 

Third,  the  maximum  growth  rate  for  each  zone  occurs  at  a value  of  q 

corresponding  to  (b  +f,+)jla)  [Figs-  6(a)  and  (b)  ] . Finally,  when 

9 

instability  does  exist,  the  maximum  growth  rate  can  be  a substantial 
fraction  of  ion  cyclotron  frequency 
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Shown  in  Fig.  7(a)  for  1=2  and  in  Fig.  7(b)  for  l=b  are  plots  of 
normalized  growth  rate  Y/“ci  versus  the  reversal  parameter  n [Eq.  (97)] 
fora=0.2  and  parameters  otherwise  identical  to  Fig.  6(b).  Comparing 
Figs.  7(a)  and  6(b),  we  note  that  the  stability  properties  are  almost 
identical  for  ct=0  and  o=0.2,  although  the  maximum  growth  rate  does 
decrease  slowly  as  ct  increases.  We  therefore  conclude  that  the 
influence  of  plasma  outside  the  layer  on  stability  behavior 
is  weak,  at  least  when  a is  sufficiently  small  (a<0.5,  say). 

Moreover,  we  also  note  from  Fig.  7 that  the  width  of 
the  instability  zones  is  reduced  for  increasing  values  of  azimuthal 
harmonic  number  l.  Furthermore,  high-harmonic  perturbations  are  easily 
stabilized  as  q approaches  zero. 

The  dependence  of  stability  properties  on  plasma  density  is 
illustrated  in  Fig.  8(a)  for  1=2  and  in  Fig.  8(b)  for  £=4,  where  the 
normalized  growth  rate  is  plotted  versus  Jl^/c  for  rt=0.82 

and  parameters  otherwise  identical  to  Fig.  7.  Note  that  the  growth 
rate  is  a decreasing  function  of  m^Rg/c.  This  feature  is  also  evident 
from  Eqs.  (87)  and  (90)  for  a=0.  We  further  note  from  Fig.  8(b) 
that  the  system  is  completely  stabilized  above  some  critical  value  of 
wpiR0^C'  For  examP^e»  from  Fig.  8(b),  the  l=k  perturbation  is  stable 
for  Wp^Rq/c>7.  In  this  regard,  we  conclude  that  the  layer-plasma 
configuration  can  be  completely  stabilized  provided  the  plasma  density 
is  sufficiently  high. 

Of  considerable  interest  for  experimental  application  is  the  stability 
behavior  for  a field-reversed  configuration  with  n<0.  Figure  9 shows 
a plot  of  normalized  growth  rate  versus  n obtained  from  Eq.  (97)  for 
SL=2  [Fig.  9(a)]  and  £=4  [Fig.  9(b)],  and  equilibrium  parameters 
Wp^R^/c^lO,  a=0.2,  Rq/Rc=0.5,  and  a/Rg«0.05.  As  discussed  in  Sec.  V.B, 


for  ri<0,  Eq.  (97)  is  valid  only  when  -l<n<-2.  /(l+2i.  ; [Eq.  (82)], 


and  |k|  p_»l.  Moreover,  the  electromagnetic  coupling  parameter 

r kq 

|k|  __  decreases  to  zero  when  n approaches  minus  unity.  In  this  regard, 

r_K0 

the  plots  in  Fig.  9 are  presented  only  for  the  range  -0.85<n<-0.5. 

We  note  from  Figs.  9(a)  and  9(b)  that  the  instability  growth  rate 
decreases  considerably  as  |n|  approaches  0.5.  The  stabilization 
for  low  value  of  [n|  is  associated  with  finite  beta  electromagnetic 
effects.  Moreover,  this  stabilization  is  most  pronounced  for  high 
azimuthal  harmonic  numbers  [compare  Figs.  9(a)  and  (b) ] . We  conclude 
from  Figs.  7 and  9 that  the  system  is  most  stable  when  the  magnetic 
compression  ratio  n approaches  zero. 

Finally,  we  conclude  this  section  by  pointing  out  two  areas 
in  which  the  present  analysis  can  be  extended  in  a relatively 
straightforward  manner.  First,  the  analysis  can  be  generalized 
to  allow  for  z-dependent  perturbations.  Second,  stability  properties 
can  also  be  investigated  for  arbitrary  plasma  density  [i.e.,  arbitrary 
value  of  | K | defined  in  Eq.  (71)]  by  solving  the  eigenvalue  equation 
(64)  numerically  in  parameter  regimes  where  the  solution  is  not  accessible 
analytically.  These  two  generalizations  are  currently  under  investigation 
by  the  authors. 
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VI.  CONCLUSIONS 


In  this  paper,  we  have  investigated  the  electromagnetic  stability 
properties  of  an  intense  P-layer  immersed  in  a dense  plasma  background. 

The  equilibrium  and  stability  analysis  was  carried  out  within  the 
framework  of  a hybrid  Vlasov-fluid  model  in  which  the  background  plasma 
electrons  and  ions  are  described  as  macroscopic,  cold  fluids  confined 
by  the  axial  magnetic  field  B^(r)^z,  and  the  layer  ions  are  described 
by  the  Vlasov  equation.  Moreover,  the  equilibrium  and  stability  properties 
were  calculated  for  the  case  in  which  the  background  plasma  has  a 
step-function  density  profile  and  the  layer  ions  are  described  by  the 
rigid-rotor  distribution  function  in  Eq.  (2).  Various  equilibrium 
properties  were  calculated  in  Sec.  II.  One  of  the  most  important 
features  in  the  equilibrium  analysis  for  a thin  layer  is  that 
Budker’s  parameter  v for  the  layer  ions  is  directly  related  to  the  magnetic 
compression  ratio  n by  v=(l-n) /(1+n) . Electromagnetic  stability 
properties  were  investigated  in  Secs.  III-V,  assuming  that  all 
perturbed  quantities  are  independent  of  axial  coordinate  (3/9z=0). 

A formal  stability  analysis  was  carried  out  in  Sec.  IV.  Equation  (69), 
when  combined  with  Eq.  (68),  constitute  one  of  the  main  results  of 
this  paper  and  can  be  used  to  investigate  stability  properties  for 
a broad  range  of  system  parameters.  A detailed  analytic  and  numerical 
investigation  of  electromagnetic  stability  properties  was  carried  out 
in  Sec.  V for  a dense  background  plasma.  It  was  found  that  the 
effects  of  (a)  the  equilibrium  magnetic  field  depression  produced  by 
the  P-layer,  (b)  transverse  magnetic  perturbations  (6B#0) , (c)  small 
(but  finite)  transverse  temperature  of  the  layer  ions,  and  (d)  the 
dielectric  properties  of  the  background  plasma,  all  have  an  important 
influence  on  stability  behavior.  For  example,  for  a dense  background 


1 

plasma,  the  system  can  be  easily  stabilized  by  a sufficiently  large 
transverse  temperature  of  the  layer  ions . 
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FIGURE  CAPTIONS 


Fig.  1 
Fig.  2 
Fig.  3 
Fig.  4 


Equilibrium  configuration  and  coordinate  system. 

Density  profile  [Eq.  (9)]  for  the  layer  ions. 

Axial  magnetic  field  profile  [Eq.  (14)]. 

Sketch  of  the  envelope  function  ij>(r)  versus  r [Eq.  (11)] 


with  i(i(R  )=0=i];(R„)  and  [3i|>/9r]  _D  =0  . 
1 1 r-R0 


Fig.  5 


Electron  and  ion  density  profiles  [Eq.  (44)]  for  the  background 
plasma. 


Fig.  6 


(a)  Plot  of  sum  of  wave  admittances  (b_+b+)^_£^  [Eq.  (96)] 

0 


versus  n,  for  1=2,  a>  ,Rri/c=4  and  a=0.  (b)  Plot  of  normalized 

pi  u 


growth  rate  y/u  . versus  n [Eq.  (97)]  for  R./R  =0.5,  a/R„=0.05, 
Cl  u c u 


and  parameters  otherwise  identical  to  Fig.  6(a). 


Fig.  7 


(a)  Plot  of  normalized  growth  rate  y /£>ci  versus  n [Eq.  (97)] 


for  SL=2,  a=0.2,  and  parameters  otherwise  identical  to  Fig.  6(b) 


(b)  Plot  of  normalized  growth  rate  y/ui^  versus  n [Eq.  (97)] 


for  1= 4 and  parameters  otherwise  identical  to  Fig.  7(a). 


Fig.  8 


Plots  of  normalized  growth  rate  y/ui  . versus  w .Rn/c  [Eq.  (97)] 

ci  pi  u 


for  (a)  1=2  and  (b)  1= 4,  with  q=0.82  and  parameters  otherwise 
identical  to  Fig.  7. 
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